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The aim of this study was to determine the strengths of the couphng of electrons with the polar long- 
wavelength transverse optical (TO) vibrations from infrared spectroscopy data. This determination 
is made by means of a simple relationship between the electron-TO-phonon interaction constant 
and material parameters, based on a parametrization of the electron-TO-phonon coupling in terms 
of the long-range dipole-dipole interaction. The combination of experimental data employed here 
allowed us to calculate directly the relevant constants for a number of representative polar insulators 
, and to show that in ferroelectrics the interband electron-TO-phonon interaction at the F point is 

04 ■ essentially strong. In these calculations, infrared spectroscopy methods proved to be an effective 

(— I ' tool for study of the main properties of electron interaction with polar long-wavelength TO phonons. 
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I. INTRODUCTION 



In solid state science, there are a number of techniques for probing detailed characteristics related to the structural 
and dynamical properties of materials, which are therefore capable of numerical prediction of electron-phonon coupling 
constants from experimental data. In particular, methods of vibrational spectroscopy have been shown to yield reliable 
estimates of several fundamental microscopic parameters, including relevant constants of electron-phonon interactions 
(e.g., flUl'l and references therein). 
^ Two common properties of the lattice dynamics of a polar crystal are (i) the existence of dipole moments associated 
r"| with optical vibrational modes and (ii) the difference in long-range fields given by longitudinal and transverse optical 
O , modes at the zone center '6^. As shown in our previous study [vj, in polar materials, the coupling of electrons with 
O the polar long-wavelength transverse optical (TO) vibrations is essentially influenced by the features of the Coulomb 
interaction between electrons and the lattice ions. In turn, these features are directly related to the presence of 
the dipole moments of the TO vibrational modes. At the same time, dependencies of the dipole moments on the 
I> ' ionic degrees of freedom give rise to infrared (IR) spectra [1, Q , which makes IR spectroscopy an important tool for 
CO investigating the low- frequency dynamics of the IR-active TO modes in polar materials (e.g., Il0]). 

■ The characteristics of the IR spectra can be understood in terms of effective charges belonging to the TO vibrational 
I modes. The relative values of these charges are directly related to the strength of the corresponding long-range 
electron-TO-phonon couplings '?']. By combining these features, we see that information on the intensities of IR- 



active vibrations provided by IR spectroscopy measurements can be successfully used for quantitative modeling of 



the interactions of electrons with zone-center TO phonon modes. Thus, the analysis of the IR spectra data in terms 
of the dielectric function can serve as an alternative to ab-initio methods to obtain from an experiment independent 
estimates of the electron-TO-phonon (el-TO-ph) coupling constants in a given material. 

To the best of our knowledge, with only a few exceptions, no direct estimates of the el-TO-ph interaction constants 
were previously available to provide numerical values for most polar materials, especially for ferroelectric systems in 
which the strength of the el-TO-ph interaction should be especially strong In particular, even for several well- 
ed ' studied compounds, such as BaTiOa, SrTiOa, and SnTe, there was little information on how different or similar the 
relevant constants of the electron-TO phonon coupling are. In this work, we have performed a numerical investigation 
of the el-TO-ph interaction constants for a number of representative polar compounds using the theoretical scheme 
proposed in [tJI. To obtain a more detailed picture, we present a comparative analysis of the results of our calculations. 
As expected [7| j very large values of the interband el-TO-ph interaction constants were found to be a general feature 
of ferroelectric materials. 
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II. THEORETICAL DETAILS 



Contributions of the IR active optical phonons to the dielectric function e{uj) of the far-IR spectral range, which 
are the subject of IR spectroscopy measurements and analysis, are usually modeled by a classical damped harmonic 
oscillator model El : 

e{iu) = eoo + 5^ ^2 "^2^^ ■ (1) 

j 3 'J 

with rij, 7j and S{j) the j-th zone-centre TO vibrational mode frequency, damping and dipole oscillator strength, 
respectively; eoo is the core (electronic high-frequency) contribution to the dielectric function. 

In our previous work 0, within a first- principles methodology, we predicted theoretically the significant increase 
of the constants of the el-TO-ph coupling in ferroelectric materials, showed how the el-TO-ph interaction constants 
might be dependent on macroscopic material parameters, and obtained analytical equations allowing us to estimate 
el-TO-ph interaction strengths in a wide range of polar dielectrics. In the limit of a zero phonon wave- vector, the el- 
TO-ph interaction represents the long-range electron-lattice dynamic hybridization of the electronic bands of opposite 
parities, which takes into account the relevant s-, p- and d-channels (e.g., p^ - fisj and references therein): 



H^i-vh = N-'^' E E SL' C^) «^k«.'k yo, . (2) 



Here the matrix elements of the cl-TO-ph interaction at the F point, g^^, (k) , serve as the central parameters of the 
effective Hamiltonian ([2]), which characterize the interband scatterings of the bond electrons due to the zone-center 
TO phonons; a+ (a) are the creation (annihilation) operators for electronic states | ak > of the valence and conduction 
bands (u, a ) in an insulator, yoj denotes the zero-momentum TO mode normal coordinates, and N is the number 
of unit cells. For instance, as applied to a family of high polar crystals such as the ferroelectric ABO3 perovskite 
oxides, hybridization ([2|) involves significant mixing between the filled O 2p and the empty d° (Tif+, Nb^+, Zr''+, 
Ta^+, Mo^+, W^+, etc.) electronic states caused by the IR-active TO Fi„ soft vibrations [11, 

For the last several decades, the Frohlich-type mod el (f2| was a useful prototype for intensive studies of various 
properties of both typical insulating perovskites [l3l.[l6l [T8l - i21j and the A^^B^^ narrow-gap semiconductors and their 
alloys (23 - |26j . It was recently demonstrated to be applicable to give a theoretical description of the ferroelectricity 
found in BiFeOa-type multiferroics [13, . In light of these studies and of the fundamental significance of electron- 
phonon coupling in the understanding of various physical properties of polar compounds (e.g., ,29i] ) . the essential 
question then naturally arises of what method could be applied to estimate (regardless of a model used) the numerical 
values of the el-TO-ph interaction constants gj in a polar material. As shown in Ref. [3] , one particular advantage 
in the study of TO lattice vibrations is that one can associate the interaction between electrons and the polar long- 
wavelength TO phonons with the long-range dipole-dipole interaction. This link, which was not anticipated by 
previous theoretical considerations, provides the following relationship 

So, = ^E,M,S{j), ^ (3) 

Here the first equation relates, for each zone-center TO vibration of the branch j and of the reduced mass Mj, the 
bare constant of the el-TO-ph coupling at the F point, goj, with the macroscopic material constants, such as S(j), 
£00 , and a bond-gap energy Eg ; gQj represents the effective k- independent interband el-TO-ph interaction defined by 
a proper average of the squared el-TO-ph matrix elements over the given electronic states: 

^0. = E ^E I sL' (k)i' \EA^)'-E.m ' 

(T > (7 k 

where Ea(k) are the energies of the Bloch states of wavevector k in the crth band. The second equation in defines 
the renormalized el-TO-ph coupling constant gj that accounts for local-fields effects-induced partial screening of the 
bare el-TO-ph interaction. 

With the help of Eqs. ([3]), one can directly evaluate the el-TO-ph coupling constants via the dipole oscillator 
strengths. The dipole oscillator strengths can in turn be readily determined through IR spectroscopy data analysis. 
Therefore, employing the relevant information from the IR spectra concerning the IR-active optical phonons and the 
dielectric function behavior in the far-IR spectral range, and using Eqs. ([3]) in combination with the experimental 
data for eoo and Eg, enables us to determine numerical values of the el-TO-ph coupling constants for polar systems 
of interest. 
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TABLE I: The calculated values of the el-TO-ph interaction constants provided for each polar TO mode of a cubic phase of the 
ferroelectric ABO3 perovskite oxides and A^^B^^ narrow-gap semiconductors. The material parameters used are also given. 
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III. RESULTS AND DISCUSSION 

We performed numerical calculations with Eqs. ([3|) for representative polar compounds belonging to several classes 
of materials: typical insulating perovskites, the A^^B^^ narrow-gap semiconductors, and various non-ferroelectric 
wide-gap binary compounds. The evaluation of goj and gj for these compounds involved the knowledge of material 
parameters (<S'(j), £00, Eg) whose experimental values were taken from the literature or, when necessary, calculated 
from the corresponding experimental data. Lists of all compounds considered and the calculated el-TO-ph coupling 
constants and material parameters used are given in Tables 1 and 2. 

Additionally, in both tables we presented the values of the dimensionless coupling parameter Kj — gj /{EgKo) cal- 
culated for low-frequency TO modes. Here Kq — Ane^/vo, and Vo is the cell volume with respect to the formula 
unit. From a physical point of view, Kj is represented as a ratio of the electronic contribution g^/Eg, determined by 
the electronic structure and the electron-ion potential, and some lattice force constant Kq. This specially introduced 
parameter is of considerable interest due to its direct ability to supply a measure of the zone-center TO phonon 
softening arising from the interband scatterings of the bond electrons. As seen from Tables 1 and 2, only the wide-gap 
ferroelectric materials are characterized by the high values of the parameter k. Since in the cubic binary crystal 
systems one can set j = 1, we also omitted the indices j throughout Table 2. Note that for the most part of binary 
compounds represented in this table, we calculated the dipole oscillator strengths S with the following correspon- 
dence: i^j^o ~ ^To) ~ ^ which relates the squared frequencies of the longitudinal and transverse zone-center 
optical vibrations, fi^Q and ^xo^ respectively, with S. 

Based on the comparative analysis of the results obtained and summarized in both tables, we can draw several 
general conclusions and discuss the role of the el-TO-ph coupling related to the ferroelectricity phenomenon: 

i) A simple calculation recently proposed in [Tj as a quantitative measure of the el-TO-ph interaction is shown to be 
applicable for a wide number of selected polar insulators and semiconductors, including representative ferroelectric 
compounds. 

ii) Among the compounds considered, it was clearly demonstrated that the el-TO-ph coupling strength is small in alkali 
halides, intermediate in binary oxides, both the alkaline earth and 3d-transition metals, and, in contrast to the first 
two cases, significantly large in ferroelectric perovskites. Notice also that the values of the el-TO-ph coupling constants 
gj tabulated in Table 1 for A^^B^^ compounds agree with the relevant values from semi-empirical estimates [l2j. The 
values of gj given in Table 1 for the wide-gap ferroelectrics can be compared with the only known semi-empirical 
estimate of Ref. [l2| for BaTiOa , which, in turn, should be recast by applying the accurate magnitude of the effective 
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TABLE II: The calculated values of the el-TO-ph interaction constants in a number of cubic binary crystal systems. The 
material parameters used are also given. 
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charge: using the value 8.33 estimated from the data of Table 1, one obtains the corrected estimate 4.2 eV/A, which 
is close to the value 4.7 eV/A of Table 1 for BaTiOa. 

iii) There is a nontrivial one-to-one correspondence with the zone-center TO vibrational mode effective charge Z*, 
or equally with the Born effective charge. It is seen that the charge Z* sets the scale of the el-TO-ph coupling; for 
instance, in pure ionic wide-gap materials, where the Born effective charges are very close to the nominal rigid-ion 
values, the actual el-TO-ph coupling strength is small. To further quantify the interplay between these physical 
quantities, we plot their interrelation in Fig. 1, which uniquely displays the substantial growth of the el-TO-ph 
coupling with a rise in the effective charge: k cx |^*|^, jus t as predicted Due to the anomalously large Born 
effective charges, which occur in ferroelectric compounds |30l. Is^ . a marked increase in interband hybridization leads 
to a drastic enhancement of the interband el-TO-ph scatterings. In the limit of pure ionic wide-gap compounds, 
multiband channels of electronic scatterings associated with zone-center TO phonons turn out to be marginal, as the 
electronic polarizability effects are not high enough to exhibit a noticeable el-TO-ph hybridization between different 
bands. Thus, our calculations show that strong el-TO-ph coupling should be a general feature of wide-gap compounds 
with high electronic polarizabilities. 

As a matter of fact, the same results hold in the case of A^^B^^ narrow-gap semiconductor systems where, however, 
one should take into account two factors: smallness of Eg and sufficiently high range of eoo ■ It is easy to see from 
Table 1 that the relevant values of the product MjS{j) for both wide-gap and narrow-gap ferroelectrics appears to 
be close to each other. At the same time, in non- ferroelectrics, as seen from Table 2, these values are significantly 
smaller. Moreover, for both groups of ferroelectrics, the unscreened dimensionless el-TO-ph coupling parameter, 
Kjicoo + 2)/3, as seen from Table 3, turns out to be about the same order of magnitude. Meanwhile, for the non- 
ferroelectric compounds considered, the values of this parameter differ significantly and are much smaller than in the 
ferroelectrics. 
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TABLE III: The calculated values of the unscreened dimensionless el-TO-ph coupling parameter for a selected number of polar 
compounds in a cubic phase. 



KNbOs 


BaTiOa 


PbTe 


SnTe 


PbSe 


PbS 


BaO 


MnO 


SrO 


TlCl 


CsCl 


LiF 


Kj (^-+2' 7.3 


5.8 


5.8 


4.3 


2.4 


1.7 


0.61 


0.53 


0.46 


0.36 


0.14 


0.08 



IV. CONCLUSION 

To summarize the results from this study, we refer to our previous work [3| , where it was shown that the zone-center 
TO vibrational mode effective charge Z* can be considered to be the key macroscopic parameter of the el-TO-ph 
coupling strength. As is well-known, the Born effective charges play a fundamental role in determining the lattice- 
dynamical properties of insulating crystals and are a powerful tool for investigating the dielectric and ferroelectric 
properties of materials. In particular, they measure the coupling of IR radiation to the optic phonons. In materials 
where the el-TO-ph coupling is operative, it can be verified by spectroscopy measurements of the IR-active TO mode. 

In this paper, we pioneered the application of experimental data to estimate the interband el-TO-ph interaction 
constants in a number of representative polar compounds. We sustained theoretical arguments [7!] with numerical 
calculations of the interband el-TO-ph interaction constants and the further comparative analysis that the large 
interband el-TO-ph interaction is a special microscopic feature of the ferroelectric materials. Thus, one can conclude 
that the presence of anomalously large Borns effective charges is the key signature that the interband el-TO-ph 
coupling is essentially strong in a ferroelectric material. In contrast, in non-ferroelectrics, the strength of the el-TO- 
ph interaction is not necessarily high enough due to their lower polar nature. This result corresponds to the main 
assumption of the vibronic theory [la - llSt regarding the existence of sufficiently strong interband el-TO-ph coupling in 
displacive ferroelectrics. The strength of the interband el-TO-ph interaction can therefore serve as a direct indicator 
of the extent to which a crystal lattice is close to a possible ferroelectric instability. 
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Figure captions 

Fig. 1: Dependence of the square root of the dimensionless el-TO-ph couphng parameter on zone-centre TO 
vibrational mode effective charge. Triangles shown on the figure are the values of \Z*\ and i/k determined separately 
for each compound. The values of \Z* \ are calculated from the experimental data given in Tables 1 and 2. The direct 
line is the resulting fit. 
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